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We determine the relative branching fractions of semileptonic B decays to charmed final states. 
The measurement is performed on the recoil from a fully reconstructed B meson in a sample of 362 
million BB pairs collected at the T(4S) resonance with the BABAR detector. A simultaneous fit to 
a set of discriminating variables is performed on a sample of B — ► DXl~vt decays to determine 
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the contributions from the different channels. We measure T(B — * Di ve)/T(B — ♦ DXl v£) = 
0.227 ± 0.014 ± 0.016, T(B~ -» D*£~u e )/T(B~ -> DXi'vi) = 0.582 ± 0.018 ± 0.030 and F(B~ -c 
D**£-v e )/r(B- -» DXl'ue) = 0.191 ± 0.013 ± 0.019 for the charged B sample, and r(B« -» 
D£-vi)/r(W -f DXl'Dt) = 0.215 ±0.016 ±0.013, r(B° -> D*£~i?t)/T(B° -> DXl~v t ) = 0.537± 
0.031 ± 0.036 and r(B° -» D** T vt) IT(FP> -> DXl~v t ) = 0.248 ± 0.032 ± 0.030 for the neutral B 
sample, where uncertainties are statistical and systematic, respectively. 



PACS numbers: 13.20He,12.38.Qk,14.40Nd 

The determination of exclusive branching fractions of 
B — > X c l~vi decays is an essential part of the B-factory 
program to understand the dynamics of 6-quark semilcp- 
tonic decays and to determine the relevant Cabibbo- 
Kobayashi-Maskawa matrix elements [l|. The mass of 
the hadronic system X c , recoiling against the leptonic 
pair, is a crucial observable both in the extraction of 
\V c i,\, in exclusive semileptonic decays, and in isolating 
B — > X u l~vt decays to determine \V u b\- It is also needed 
for the measurement of heavy quark masses and other 
non-perturbative OPE (Operator Product Expansion) 
parameters from the distribution of spectral moments. 
This mass spectrum can be better understood by a study 
of the yields of the different D meson states in semilep- 
tonic decays. Current measurements [1, [H, 0, [|[ show a 
possible discrepancy between the sum of exclusive rates 
and the inclusive semileptonic decay width Q. While 
B — ► D£~p£ and B — > D*l~vg decays account for about 
70% of this total, the contribution of other states, in- 
cluding resonant and non-resonant D^tt decays, is not 
yet well measured and is a possible explanation of this 
discrepancy. 

In this paper, we present a novel technique to extract 
the exclusive relative branching fractions for B — > D£~Dg, 
B -> D*l~vi and B -^_D**e~Pt, with I = e, fx Q, from 
an inclusive sample of B — > DXi~Dg events, where X can 
be either nothing or any particle (s) from a semileptonic B 
decay into a higher mass charm state, or a non-resonant 
state. We denote by D** any hadronic final state, con- 
taining a charm meson, with total mass above that of the 
D* state, thereby including both Dj excited mesons and 
D^+nir non-resonant states. This technique ensures 
sensitivity to all hadronic final states containing a D me- 
son, thus helping us to understand the role of excited D 
states in saturating the inclusive semileptonic rate. 

This analysis is based on data collected with the BABAR 
detector Q at the PEP-II asymmetric-energy e + e~ stor- 
age rings. The data correspond to an integrated lumi- 
nosity of 339.4 ftr 1 recorded at the T(4S) resonance, 
or, equivalently, about 362 million BB pairs. A detailed 
GEANT4-based Monte Carlo (MC) simulation [| of BB 
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and continuum e + e~ — * // (/ = u,d, s,c,t) events 
has been used to study the detector response and its 
acceptance. The simulation models B — ► D*i~vi de- 
cays using HQET-based calculations [lOj, B — ► Dt~vi 
and B -> D**(~>_ D^w)£~v e decays using the ISGW2 
model [HJ, and B — > D^irl~Dg decays using the Goity- 
Roberts model [HI]. 

We select signal _B-meson decays in events containing 
a fully reconstructed B meson (B tag ), which allows us 
to constrain the kinematics, to reduce the combinato- 
rial background and to determine the charge and fla- 
vor of the signal B. We choose a set of three largely 
uncorrelated variables to discriminate between the dif- 
ferent semileptonic decay modes in the reconstructed 
B — > DX£~i>£ sample. These are: i) the lepton mo- 
mentum in the center-of-mass (CM) frame, \pi\] ii) the 
missing mass squared reconstructed with respect to the 
Di system, which corresponds to the mass of the Xvi sys- 
tem, m 2 miss D = (p r -pB tag ~Pd -Pe) 2 , where p { is the 
four momentum in the CM frame of the reconstructed 
state i; and iii) the number of reconstructed charged 
tracks in addition to those used for reconstructing the 
Di system and the B tag , N tr ks- In order to reduce the 
sensitivity to the modeling of the decays to the different 
charm states, the shapes of these variables are extracted 
from data, using exclusive samples highly enriched in the 
relevant decay modes. The relative D, D* and D** con- 
tributions are then determined by a multiparameter fit 
to the inclusive sample. 

We select semileptonic B decays that contain one fully 
reconstructed D meson and that recoil against a fully 
reconstructed B tag decaying hadronically. To obtain a 
high reconstruction efficiency, the analysis exploits the 
presence of two charmed mesons in the final state: one 
used for the exclusive reconstruction of the B tag , and 
another in the semileptonic B decay. 

The event reconstruction starts from the semileptonic 
B decay, selecting a charm meson and a lepton with mo- 
mentum in the CM frame higher than 0.6 GeV/c and the 
correct charge-flavor correlation. Candidate D° mesons 
are reconstructed in the K~ir + , K^k+it , K~ir + it + ir~ , 
K a s TT+TT-, ifgTT+TT-Tr , K° s tt , K+K~, tt+tt", and K° S K% 
channels, and D + mesons in the K~tt + it + , K~ir + ir + ir° , 
A'°tt+, K° s n + n , K+R-1T+, K° S K+, and A°7r+7r+7r- 
channels. In events with multiple candidates, the candi- 
date with the largest D-i vertex fit probability is selected. 
We then select a fully reconstructed B tag meson candi- 
date. We reconstruct B tag decays of the type B — > DY, 
where Y represents a collection of hadrons with a total 
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charge of ±1, composed of rii7r ± + n 2 K ± + n 3 Kg + n 4 7r°, 
where n\ + n 2 < 5, n 3 < 2, and n 4 < 2. Using D°(D + ) 
and Z?*°(L>* + ) as seeds for B~(B a ) decays, we recon- 
struct about 1000 different decay chains. 

The kinematic consistency of a B tag candidate with a 
JS-meson decay is checked using two variables: t he beam- 
energy substituted mass rriEs = -\/s/4 — |ps| 2 , and the 
energy difference AE = EB — yfs/2- Here y/s refers to the 
total CM energy, and \ps\ and Eb denote the momen- 
tum and energy of the B tag candidate in the CM frame. 
For correctly identified B tag decays, the tties distribution 
peaks at the B meson mass, while AE is consistent with 
zero. We select the B tag candidate that has no daughter 
particles in common with the charm meson and the lep- 
ton from the semileptonic B decay, m es within the signal 
region defined as 5.27 GeV/c 2 < m ES < 5.29 GeV/c^, 
and the smallest |A.E| value. Mixing effects in the B° 
sample are accounted for as described in [l3j |. 

The B — ► DXi~Di decays are identified by relatively 
loose selection criteria. We require the reconstructed 
ground-state charm meson invariant mass M D o (M D +) 
to be in the range from 1.850 (1.853) GeV/c 2 to 1.880 
(1.883) GeV/c 2 and the cosine of the angle between the 
directions of the D candidate and the lepton in the CM 
frame to be less than zero, to reduce background from 
non-B semileptonic decays. 

After these selection criteria, the sample contains lep- 
tons from prompt B decays, as well as cascade B de- 
cays, in which the lepton does not come directly from 
the B. There are also background sources of leptons, 
such as photon conversions and Dalitz 7r° decays, com- 
binatorial BB background and continuum events, that 
need to be subtracted. The contamination from cascade 
B decays, about 15.1 (17.8)% of the total B~(B°) sam- 
ple, is subtracted using the simulated MC distributions 
for these backgrounds. These events are reweightcd to ac- 
count for differences among the branching fractions used 
in our MC simulation and the latest experimental mea- 
surements The photon conversion and n° Dalitz 
decay backgrounds (less than 0.8% of the total electron 
sample) are removed using a dedicated algorithm, which 
performs the reconstruction of vertices between tracks of 
opposite charges whose invariant mass is compatible with 
a photon conversion or a ir° Dalitz decay. The contribu- 
tions of combinatorial and continuum B tag backgrounds 
are estimated from the ities sideband region 5.21 GeV/c 2 
< mss < 5.26 GeV/c 2 . The tties distribution is fitted 
by the sum of a Gaussian function joined to an expo- 
nential tail [HI for the signal and an empirical phase- 
space threshold function [l6j for the background. Cross- 
feed effects, i.e. Bf (B° ta g) candidates erroneously re- 
constructed as a neutral (charged) B, are corrected us- 
ing MC simulations. We estimate the fraction of cross- 
feed events in the reconstructed B~(B°) sample to be 
6.8% (8.1%). A total of 6396±251 (2981 ±122) events are 
selected, with an estimated purity in B~(B°) — > DX£~vi 
of 72% (73.8%). 

Exclusive samples enriched in Dt~vi, D*£~Di and 



TABLE I: Invariant mass ranges for D*° and D* + selection. 

Mode Selection Criteria 

D*° -> D°ir° 0.139 < M(D*°) - M(D°) < 0.145 GeV/c 2 

D*° -> D°7 0.133 < M(D*°) - M(D°) < 0.151 GeV/c 2 

D* + -> D°ty+ 0.141 < M(D*+) - M(D°) < 0.149 GeV/c 2 

D*+ -> D + tt° 0.138 < M(D*+) - M(D+) < 0.143 GeV/c 2 



~vi are then selected. Contributions from other 
semileptonic B decays into charm final states, where 
one or more particles from a higher mass charm state 
are missing (feed-down) or random particles are erro- 
neously associated with the charm candidate (feed-up) 
are removed. This is done by selecting signal regions in 
the missing mass squared distributions m 2 . „, = 

(pr — PBtag — ?>£>(•.**) ~ Pi) 2 corresponding to the exclu- 
sive decay being reconstructed. The selection criteria are 
chosen to maximize the sample purity. We select D* + 
and D*° candidates by requiring the invariant mass dif- 
ference between the D* and the D to satisfy the selection 
criteria in Table fl] 

The B~ -> D*°£~V£ and W -> D*+£-p e decays are se- 
lected by requiring the missing mass squared TnP,. D , a 
to be between -0.35 GeV 2 /c 4 and 0.5 GeV 2 /c 4 'and 
l m miss d*+\ to be smaller than 0.55 GeV 2 /c 4 , respec- 
tively. Feed-down events from decays to D** states are 
removed by requiring m 2 misSt D „ a and m 2 mis ^ D „ + to 

be_ incompatible with zero. The B~ — > D°l~Vi and 
B° — > D + l~T>i decays are selected by removing feed- 
down events from D* and D** states. Similar selection 
criteria are applied for B~ — * D**°£~Dg decays, with 
D **o ^ j o(*)+ 7r - j an d ~B~o _> D**+£-p e decays, with 
D**+ -> D(*)°7T+. 

The probability density functions (PDFs) of the dis- 
criminating variables, \pi\, D and N tr ks are deter- 
mined using the exclusive samples. In order to test for 
possible selection biases in the PDF shapes, the inclusive 
distributions for MC samples of B — > D£~Dg, D*l~vi and 
D^7r£~De events have been compared to those obtained 
after the exclusive event selection. Good agreement is 
found after accounting for the residual background from 
feed-down and feed-up from other modes. The PDFs 
are parameterized as sums of analytic functions, such as 
Gaussians and polynomials, with the exception of Ntrks 
which is described using histograms. 

The relative fractions of D, D* and D** decays in the 
selected inclusive sample of B — > DX£~Di events are 
obtained by a simultaneous \ 2 n ^ to the inclusive and 
exclusive \pe\, D and Ntrks distributions. The rel- 

ative fractions are floated, constraining their sum to be 
one, together with the parameters of the functions de- 
scribing the shapes of the discriminating variables. This 
results in a 35-parameter fit, which ensures that statisti- 
cal correlations between the different samples are prop- 
erly taken into account and the uncertainties in the ex- 
clusive shapes, obtained from samples of significantly 
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Lepton Momentum (GeV/c) Missing Mass Squared (GeV 2 /c 4 ) Additional Charged Track Multiplicity 



FIG. 1: Fitted \pe\ (a,b), m^ niss D , (c,d), and N tr ks (e,f) distributions for B° — > DX£ vi (top) and B — > DX£ vg (bottom). 
The PDFs corresponding to the different exclusive components are superposed with different filling styles. 



smaller size compared to that of the inclusive sample, 
are correctly propagated into the statistical uncertain- 
ties on the D, D* and D** relative fractions. Since this 
analysis does not reconstruct D** states with neutral 
pions, the N tr k s distribution for states with the same 
charged-track multiplicity is used to model these decays: 
e.g. the B~ — > D*°£~Di N tr k s distribution is used for 
modeling D (—* D*°tt°)£~ decays. For the modes 
involving a soft charged pion, such as B° — > D* + £~9i, 
the MC prediction for the additional charged-track mul- 
tiplicity distribution is used to account for inefficien- 
cies in the reconstruction of the low-momentum particle. 
MC studies show that the PDFs for the B -> D**£~v t 
component, obtained by the exclusive reconstruction of 
B — > D^ir£~Di decays, can also be used to parameterize 
B — ► D^mr£~i>i decays in the inclusive B — > DX£~i>t 
sample. The fit also accounts for feed-down and feed-up 
decays in the exclusive shapes, fixing the relative con- 
tributions to the predictions from the simulation. The 
fit performance has been extensively tested using simu- 
lated samples with varying fractions of the different decay 
modes. These tests show that the procedure adopted in 
this analysis is able to extract the decay fractions without 
any significant bias. The statistical uncertainty obtained 
by the fit reproduces the scatter of the results from inde- 
pendent samples, where the bin contents of the distribu- 
tions have been fluctuated according to their statistical 
uncertainty. The fit results for the B° — > DX£~vz and 
B~ — > DX£~Di distributions of the three variables \pe\, 
m miss d aim N tr ks are shown in Fig. 1. The fit has a \ 2 
value of 200 for 212 degrees of freedom for the B~ sample 
and 204 for 168 degrees of freedom for the B° sample. 

Several stability checks have been performed. First 
the sample has been split into sub-samples based on the 
lepton flavor and the run period and the fit has been 



repeated for each one of them. Results are consistent 
within the statistical uncertainties. As another check, 
the B — > D£~p£ and B — > D*£~vi branching fractions 
have been determined by a binned likelihood fit to the 
m miss d ancl m rniss d* distributions respectively, where 
simulated events are used to model the shape of the miss- 
ing mass squared variables for the D, D* and D** exclu- 
sive decays and the combinatorial and continuum back- 
ground. The results are in good agreement with the rel- 
ative branching fractions obtained from the fit to the in- 
clusive B — > DX£~vi sample, once we normalize them 
to the total semileptonic B branching fraction. 

Different sources of systematic uncertainties have been 
investigated and are given in Table HT1 The first source is 
due to detector effects, where the size of the uncertainties 
in the detector response are determined from data con- 
trol samples. Uncertainties related to the reconstruction 
of charged tracks are determined by evaluating the fit sta- 
bility using different track selection criteria and by a MC 
study in which we vary the track multiplicity according to 
the tracking efficiency uncertainty. The systematic error 
due to the reconstruction of neutral particles is studied 
by varying the simulated calorimeter resolution and effi- 
ciency. The systematic uncertainty from lepton identifi- 
cation is estimated by varying the tagging efficiency by 
2% (3%) for electrons (muons) and the misidentification 
probability by f 5%. 

The second main source of systematic uncertainty is 
related to the selection of the inclusive sample. A major 
contribution is due to background processes, where the 
estimated systematic error is dominated by the uncer- 
tainty on the weighting factors used to subtract B cas- 
cade decays. The uncertainty in the subtraction of the 
background from the fully reconstructed B tag decays is 
evaluated from the differences in the shapes of this back- 
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TABLE II: Relative errors (%) in the determination of F(B — > 


£>(*,**) i y t )/Y{B^DXt 


Vl). 






d— / DO . ViQ~i, 

JD j JD — * _Ly£ Vg 


B~ [B° — > D*l~v t 


b-/b° -> D**rv t 


Tracking reconstruction 




Q C/l /O Of? 


1.3/0.3 


4.63/3.14 


Neutral reconstruction 




Pi /0 ^ 
U.OO/U.O 


0.39/0.31 


0.41/0.34 


Lepton identification 




3.46/3.24 


3.71/3.57 


3.51/3.3 


B - — > DJCE isg Backgrounds 




48 /0 78 


1.69/4.13 


5.11/7.19 


R i J) xfi i/n R pppiti Q"f~T*n pt~ i nn p~rripipTiP^7' 






1.53/2.6 


3.37/6.43 


B — > DX £~~ Cross-feed corrections 




23 /0 46 


0.13/0.56 


0.71/0.97 


B — *• '£~ i>t Backgrounds 




1.81/1.46 


1.04/1.41 


1.87/2.26 


B — *• \>t Feed-down and feed-up 


corrections 


1.99/1.49 


1.31/1.34 


1.84/2.01 


B — > D^'**^£~ U£ Cross feed corrections 




0.74/0.62 


0.1/0.23 


1.33/0.74 


\pi\ and m^ ss PDFs 




3.27/1.68 


1.06/1.81 


1.64/4.8 


N trks PDF 




0.38/0.9 


0.91/0.2 


3.68/0.86 


B -> D^nir£~ut 




0.9/0.73 


1.1/0.89 


0.89/0.72 


Total Syst. 




7.06/6.19 


5.18/6.71 


9.88/12.2 



ground in the sideband and in the signal region using 
MC predictions. The systematic error due to the uncer- 
tainty in the amount of flavor cross-feed is computed by 
varying its fraction by a conservative 30%. The corre- 
sponding systematic uncertainties are evaluated for the 
exclusive samples. The analysis, relying on decay classi- 
fication in an inclusive sample, is not sensitive, at first 
order, to reconstruction efficiencies. There remains an 
uncertainty arising from possible differences in efficien- 
cies for the various channels, which is estimated from 
simulation. 

Systematic uncertainties due to the PDFs are esti- 
mated by replacing the shapes extracted from the ex- 
clusive samples with those predicted by our simulation 
and repeating the fit. Additionally the uncertainty in 
the relative D**° DW+tt - to D**° -► D« 7r° recon- 
struction efficiency is accounted for by varying the N tr ks 
distribution for the D* component. 

Systematic effects due to B — > D^*'nn£~Dg events not 
well parameterized by the B — > D**l~vt PDFs are esti- 
mated by repeating the fit with an additional component 
for these events. The corresponding PDFs are built from 
a sample of simulated B — ► D^tt-k^vi events. The ob- 
served difference in the fit results is taken as an additional 
systematic error. 



TABLE III: Fitted ratios of branching fractions with statisti- 
cal and systematic uncertainties. 





Ratio 




B- (%) 


B° (%) 


r(S 
r(B- 


^DXt 


Vt) 


22.7 ± 1.4 ± 1.6 


21.5 ± 1.6 ± 1.3 


T(B- 

r(B- 


^DXt 


~vt) 


58.2 ± 1.8 ± 3.0 


53.7 ± 3.1 ± 3.6 


r(B- 
r(B- 


^Dxe 




19.1 ± 1.3 ± 1.9 


24.8 ± 3.2 ± 3.0 



In summary, the relative branching fractions for the 
B~ -> L>°, D*°, D**°£-9 e and B° -> D+, D*+ , 
£)**+(- y l decays have been determined by a multipa- 
rameter fit to three discriminating variables in an inclu- 
sive sample of B — > DX£~vi events recoiling against a 
fully reconstructed B meson. The results are given in 
Table IIIII Apart from possible isospin violation effects, 
which are thought to be small, these three ratios are 
expected to be equal for B~ and mesons. The re- 
sults for charged and neutral B mesons are compatible 
within their uncorrelated uncertainties. Therefore the 
relative fractions have been averaged, accounting for cor- 
related errors. The results are: T{B — > Dl~vt)jY{B — > 
DX£-ve) = 0.221 ± 0.012 (stat.) ± O.OOfiJuncorr. syst.) 
± 0.010 (corr. syst.), T(B -► D*i'i> l )/T{B -> DX£~D e ) 
= 0.572 ± 0.017 (stat.) ± 0.016 (uncorr. syst.) ± 
0.022 (corr. syst.), T(B -> D**£-pi)/T(B -► DXi~D t ) 
= 0.197 ± 0.013 (stat.) ± 0.013 (uncorr. syst.) ± 
0.012 (corr. syst.), where the first uncertainty is statisti- 
cal, the second the uncorrelated systematic and the third 
the correlated systematic error. The accuracy of these 
measurements is comparable to that of the current world 
average [6J. 
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